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We study the kinetics of diffusion-limited coalescence, 
A + A — *• A, and annihilation, A + A — > 0, in random me- 
dia consisting of disconnected domains of reaction. Examples 
include excitons fusion and annihilation in porous matrices 
and along polymer chains. We begin with an exact analysis 
of A + A — > A in a finite segment. This result is applied to 
coalescence in a random distribution of segment lengths, and 
the implications for coalescence and annihilation in percola- 
tion clusters and other confined geometries are then derived 
by means of scaling techniques. 



I. INTRODUCTION 

Diffusion-limited reactions are far less understood than 
their reaction- limited counterpart. Simple reaction mod- 
els such as diffusion-limited coalescence, A + A — > A, and 
annihilation, A + A — > 0, have attracted much recent 
attention since their kinetics in one-dimensional 

space may be analyzed exactly, and they shed light on 
other less tractable reaction schemes. Moreover, co- 
alescence and annihilation serve as models of exciton 
dynamics in several real systems. Examples include 
photoexcited solitons and polarons in MX chain com- 
pounds, such as [Pt(en)2 l[Pt(e n) 2 Cl2l(BF4)4 ((en) de- 
notes ethylenedyamine) |lj,|l6|, and fusion of photo- 
generated excitons in tetramethylammonium manganese 
trichloride (TMMC) chains |l7). Kroon and Sprik (l?]] 
review many other systems which are aptly modeled by 
diffusion-limited coalescence or annihilation in one di- 
mension. 

Most theoretical studies to date assume an infinite 
medium, however, rich kinetic behavior arises from the 
finite size of the substrate in question. For times that 
exceed the typical diffusion time needed to span the size 
of the sample, finite-size effects set in and dramatically 
affect the kinetics. We argue that when the system con- 
sists of a collection of reaction-domains of varying sizes, 
this might lead to a stretched exponential decay in the 
concentration of particles, as opposed to the usual power- 
law decay. For example, disorder in a system of polymers 
might effectively fracture the medium into separate reac- 
tion zones. Indeed, stretched exponential decay of photo- 
generated excitons in MX compound chains |15],[l6| and 
in TMMC chains j[7) is a well documented fact. 

The problem generalizes to systems of higher dimen- 
sionality. Triplet exciton annihilation in naphtalene-hs 
embedded in porous matrices of napthalene-ds has been 
studied by Kopelman Here the geometry of the re- 
action zones and their characteristic sizes might be well 



described by percolation theory. Kopelman Jl8| has stud- 
ied exciton dynamics in several systems, spanning an im- 
pressive range of different geometries. 

In this paper we deal with diffusion-limited coales- 
cence, or annihilation, in random media — media con- 
sisting of disconnected reaction domains with a given 
distribution of sizes. We begin with an exact deriva- 
tion of the kinetics of coalescence, A + A — ► A, in a 
finite segment of length L, in section II. At early times, 
t <C 1/Dcq (D is the diffusion coefficient of the parti- 
cles, and Co is their initial concentration), few reactions 
take place and the concentration is nearly constant; at 
the intermediate time regime of 1/Dcq <C t <C L 2 /D the 
concentration decays in power-law fashion, c ~ l/yDt; 
and in late times, L 2 /D <C t, the decay is exponential, 
c ~ exp[— (2tt 2 D / L 2 )t\. In section III, we use these re- 
sults to derive the kinetics of diffusion- limited coalescence 
of particles confined to a set of segments characterized by 
the distribution of their lengths, ip(L), as might be the 
case for exciton dynamics in polymer compounds. For a 
completely random distribution, ip(L) = p^ 1 cxp(— pL), 
we find that the long time asymptotic decay is stretched- 
exponential, c ~ exp(— i 1 / 3 ), similar to the long time be- 
havior in the famous problem of trapping . The earlier 
time regime is dominated by a somewhat faster strctched- 
exponential decay, c ~ exp(— t 1 ' 2 ). Coalescence and 
annihilation in percolation clusters is addressed in sec- 
tion IV. Here the situation is more complicated, for the 
reaction domains are fractal and diffusion is then anoma- 
lous p0[ ]. Our analysis consists of scaling arguments, 
motivated by the more firmly grounded results in one di- 
mension. We conclude with a summary and discussion, 
in section V. 



II. COALESCENCE IN A SEGMENT 

Consider the diffusion-limited coalescence process, A + 
A — > A, where the particles diffuse with diffusion co- 
efficient D and merge instantaneously upon encounter, 
taking place in a segment of length L. An exact anal- 
ysis is possible through the method of Empty Intervals, 
known also as the method of Inter-Particle Distribution 
Functions (IPDF) (]]|. The key parameter is E(x,t) - 
the probability that an arbitrary interval of length x is 
empty (contains no particles), at time t. E(x, t) satisfies 
the diffusion equation [p3: 



d_ 



E{x,t) 



(1) 
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The coalescence reaction imposes the boundary condition 
E(0,t)=l, (2) 

and, because E is a probability, also 

E(x,t)>0. (3) 

Finally, if there are any particles in the segment then 

E(L,t) = 0. (4) 

From E(x,t) one can compute various useful quanti- 
ties. For example, the particle concentration is 



dE{x,t) t 
C{t) = dx— [ ^ 



(5) 



and the conditional probability that the next nearest par- 
ticle to a given particle is at distance x (the IPDF), is 



d 2 

p(x,t)=c(t)- 1 —E(x,t) 



(6) 



The steady state solution of eq. (M (i.e., dE/dt = 0), 
with the boundary conditions (0), (g), and (Q), is 



E s (x) = 1 - - 



(7) 



Clearly, if there is initially at least one particle present, a 
single particle would remain at the end of the process, un- 
able to undergo any further reactions. This is confirmed 
by eq. and using (|): 



1 

c s = — 
L 



(8) 



In systems consisting of various reaction domains, this 
gives rise to a residual concentration at the end of the 
process (Appendix A). Our goal is to describe the ap- 
proach of the concentration to the residual limit. 

For the full time analysis, we rewrite the empty interval 
probability as a sum 



E(x,t) =E s (x)+ E tIsns (x, t) 



(9) 



where -Etrans symbolizes the transient behavior of inter- 
est. -Etrans satisfies the same equation as E, eq. (|l|), but 
the boundary condition (0) is replaced by 



E trans (0,i) =0 • 

A spectral decomposition then yields 



E t i£ 



E. nnx , 
a n sm(— — )e 
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where the {a n } are determined from initial conditions. 

For practical purposes, we consider the case where the 
segment is initially filled with particles at concentration 
c . We assume that the particles are randomly (Poisson) 



distributed, and hence the initial empty interval probabil- 
ity, given that the segment contains at least one particle, 
is 



and 



E(x,0) 



Etrans (^,0) = 



-c L\ 



1 — e~ c ° L 



-c L L 
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, — CqL 



(12) 
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Performing a Fourier analysis and comparing to eq. (|llj) , 
we find 



2coE 2 (cos?i7r- e CoL ) 



(14) 



" (-1 + e c ° L )mr(c 2 L 2 + n 2 n 2 ) ' 
Finally, from this and eqs. (|^) and ( |TT| ) we conclude that 
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y ( e c L _ + n 2 it 2 /c 2 L 2 ) 



(15) 



The residual concentration 1/L has been eliminated, so 
that ct(i) represents the approach to the asymptotic 
limit (in a single segment of length L, and given that 
there is at least one particle at the start). 

It is easily verified that the limit L — > oo yields the 
known result for the infinite line ]l3| : 

c oo (t) = c e t/t0 eric^/t/h , t = . (16) 

If L is finite, however, at long times, t 3> t± = L 2 /(2tt 2 D, 
the sum in ( |l5| ) is dominated by the slowest decaying 
exponential: 



c L ~ a(L)e- t/41 , a(L) 



= CoL + 1 



L(l + 7r 2 /4L 2 )e c o L -1 



(17) 



III. SYSTEMS CONSISTING OF SEGMENTS OF 
RANDOM LENGTH 

We now have all the tools for the analysis of diffusion- 
limited coalescence in a system consisting of a collection 
of segments of various lengths, taken from the distribu- 
tion ip(L). The approach to the steady state is given 
by 



c(t) = (L)- 1 c L {t)L^{L)dL 
Jo 



(18) 



where (L) = J °° Lip(L) dL is the average length of the 
segments. Here and henceforth we write c(t) instead of 
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c(t) — c(oo), for the sake of brevity. The residual con- 
centration of particles at the end of the process, c(oo), is 
analyzed in Appendix A. 

Eq. (jig), in conjunction with cq. (|l5|), gives the ex- 
act concentration decay, however, some approximations 
are necessary in order to reveal the asymptotic trends. 
Thus, we replace the exact result for cl(£) with the lim- 
iting behavior of eqs. © and (H|), but expressing the 
crossover between the two regimes in terms of the length 
A = \ / 2ir 2 Dt rather than the time t\. 



Cr(t)~i C °°® , L>>A ' 
CL[t) l a^e"*/' 1 L«A. 



(19) 



Then, 



c{t)~{L)- x {\ a{L)e- l/tl Li){L)dL 



o 

+ coo (t) ^ d&} = c< (t) + c> (t) , (20) 

where c< (t) and c> (t) represent the contribution to c(t) 
from short (£ < A) and long (L > A) segments respec- 
tively. Furthermore, a reasonable approximation to ( |l6| ) 
is: 



Coo(i) 



1 



(21) 



l/c + V27r£)< 
Essentially, at times t ^S> to = (2cqD)^ 1 the initial con 



centration has no effect and c~ 



l/V27rDi. 



A. Exponential distribution of segment lengths 

Consider first a completely random distribution of seg- 
ment lengths, tp(L) = pe~ pL . This might be the case in 
chain compounds when the concentration of disorder (de- 
fects along the polymers which discourage diffusion and 
effectively terminate the chains) is p, and the average seg- 
ment length is (L) = 1/p. We also assume that cq ^> p, 
that is, the initial average number of particles per seg- 
ment is large. The contribution from short segments can 
then be estimated using a saddle point approximation: 



:(*) 



V 1/6 e" 3Tl 



r > 1 



(22) 



where r = t/(2/TT 2 Dp 2 ) is a reduced (dimensionless) 
time. c<(£) is essentially the same as the survival proba- 
bility in the famous trapping problem |l9| , and for similar 
reasons: the stretched exponential decay results from the 
longer survival of particles in very long segments, weighed 
against the relative scarcity of such segments. 

The contribution from long segments can be computed 
analytically: 



c> (t)=c 00 (t)[l+pA(t)]e-r A W 



Then, with cq S> p, and using the approximation (|21|), 
we have 



C> (t)~V^p(l + ^)e- 2VT 



T > 1 



(24) 



The two contributions, c< and c>, become equal at 
r sa 3.5. Thus, a crossover is expected between the 
stretched exponential decay with exponent 1/2 at early 
times (but greater than 2/n 2 Dp 2 ) and a stretched expo- 
nential decay with exponent 1/3 at later times. In reality, 
the concentration decay is stretched-exponential, with a 
slowly- varying exponent; from roughly 1 atrsa 1, to 1/2 
at t w 10, to 1/3 at r > 100 (figure 1). At the very 
beginning of the process, t <C 2/n 2 Dp 2 (r <C 1), the sur- 
vival probability is practically the same as in an infinite 
system: c(t) « Coo (t) of eq. ©, or fn]). 




FIG. 1. Concentration decay in an exponential distribution 
of segment lengths. Shown is ln(c) vs. r in a log-log plot, for: 
c< (dashes), c> (dash-dotted line), and c = c< + c> (solid 
line). Slopes of 1/2 and 1/3 are shown for comparison. 



B. Power- law distribution of segment lengths 

Consider next a power-law distribution of the segment 
lengths: 



L < 6, 

( 7 -l)6T- 1 i-T L>6, 



(25) 



(23) 



where b is a microscopic cutoff length, and 7 > 1. This 
is an example of a long-tailed distribution, and it is 
also characteristic of percolation clusters at criticality. 
If 7 < 2 then the average segment length diverges (in the 
thermodynamic limit, when the size of the system is infi- 
nite) . In this case the system is dominated by the largest 
segment, which then coincides with the overall size of the 
system, L, and c(t) ~ cl(£) of Section II. 

For 7 > 2 the average segment length is finite: (L) = 
[(7 — 1) /(7 — 2)]b. The contribution from short segments 
can be written as 



3 



r A / b A 2 A A 

C< (t) = (7 - 2)P- 2 I a(-)e- u -du (26) 



Since b is a microscopic cutoff lengthscale, we assume 
that Cob < 1. In this case t ^> to = \j(2c\D) implies 

2 

A ^> b. Due to the term e~ u , the main contribution 
to the integral in ( p6f ) comes from u w 1, in which case 
a(A/u) ~ 2u/A, and the upper integration limit may be 
extended to +oo. Thus, 



(t) ~ ( 7 - 2)r(l±i, l)F- 2 (2^ J Di)( 1 -^/ 2 , (27) 



where T(a, a;) = J°° z a_1 e~ 2: is the incomplete gamma 
function 

The contribution from the large segments, for i 3> to 
(and cgfe < 1), is 



C>(t) = Coo(t) ( - 



7-2 



(28) 



Thus, the contributions from short and long segments are 
of the same order of magnitude at all times t 3> to , and 
the overall concentration decays in power-law fashion, 
c(t) - (l/b)(b 2 / Dtp- 1 ^ 2 , faster than the 1/v^-decay 
of coalescence in the infinite line. 



IV. REACTIONS IN PERCOLATION CLUSTERS 

The physical properties of three-dimensional systems 
with disorder are well represented by the percolation 
model of second-order phase transitions p2[ . A relevant 
example is triplet exciton annihilation in naphtalcnc-hs 
embedded in porous matrices of napthalene-ds jl£| . The 
domains of reaction of naphtalene-hs can then be mod- 
eled by percolation clusters. Suppose that a reaction do- 
main, i.e., a percolation cluster, has mass s. Then the 
concentration of particles in the coalescence (or annihi- 
lation) process is [|o|Jl|j2|,||: 



c s (t) 



lp-t/ti 



S > (7, 

s < a, 



(4/2 



(29) 



where d s sw 4/3 is the spectral (or fracton) dimension of 
the clusters, and a is the crossover mass between alge- 
braic and exponential decay. 

Consider first percolation well below the critical 
threshold, when the distribution of cluster masses is 



n(s) 



pe 



The typical cluster mass 1/p is smaller 



the farther one is from the percolation transition. Then, 
following a reasoning similar to the one of Section III. A, 
we find 



c(t) ~ pr 



d s ' /2(2+d s ') 



'/(2 + d B ') 



-d B /2 



r ds/2 



(30) 



where r = p 2 / ds t is the reduced time, and the first and 
second terms on the r.h.s. represent c< and c> respec- 
tively. We note that the largest (and rarest) clusters in 
percolation below criticality — which are the ones that 
actually determine the asymptotic behavior of c< — are 
different from the typical clusters. They are known as 
"lattice animals" , and their spectral dimension d s f» 1.18 
(in d = 3) is somewhat smaller than d s . Thus, the ex- 
ponent characterizing the stretched exponential decay 
of c(t) varies slowly from d s /2 ps 0.67 at r ps 1, to 
d s ' /(2 + d s ') ps 0.37 at later times. 

Exactly at criticality, the distribution of the largest 
(dominant) clusters is given by n(s) ~ s~ T , with r 
slightly larger than 2 (r ps 2.186, in three-dimensional 
percolation). Following the argument for power-law dis- 
tributions we then find 



c(t) ~ t 



-(r-l)d a /2 



(31) 



where once again the contributions from small and large 
clusters are comparable. Using the known values of the 
various exponents in three dimensions we predict an al- 
gebraic decay: c{t) ~ t~ a79 . 

Well above criticality the system is dominated by the 
spanning (infinite) cluster. This cluster has fractal struc- 
ture, and is characterized by the spectral dimension 
d s ps 4/3, only for lengthscales shorter than the correla- 
tion length £. The correlation length is infinite at critical- 
ity, but shrinks rapidly as one moves beyond the critical 
threshold. It is also the typical size of the finite clusters 
in the system. At lengthscales larger than £, the span- 
ning cluster is similar to a regular d-dimensional object. 
Thus, in three-dimensional percolation above criticality 
we find 



it) 



£-ds/2 
i- 1 



t < ^ 
* > ^ 



(32) 



where d w ps 3.88 is the walk dimension of diffusion in the 
fractal clusters. The 1/i-concentration decay at length- 
scales larger than £ takes place because the upper critical 
dimension of diffusion-limited coalescence (and annihila- 
tion) is d = 2, so that in the effectively three-dimensional 
cluster the kinetics is reaction-limited. 



V. SUMMARY AND DISCUSSION 

Starting from the exact solution of diffusion-limited co- 
alescence (A + A — ► A) in a segment of finite length L, we 
have computed the approach of the concentration to its 
steady-state limit, in a system consisting of disjoint seg- 
ments characterized by the distribution of their lengths, 
ip{L). For an exponential distribution, ip(L) = pe~ rhoL , 
as might be the case in polymer chains with an element 
of disorder, we find three different regimes: (a) an early 
time regime where the concentration decay is similar to 
that in the infinite line, c ~ 1/yDt, (b) an intermedi- 
ate time regime, dominated by a stretched-exponential 



4 



decay, c ~ e at , and (c) a late time regime where 
c ~ e~P l (a and (3 are known constants). In practise, 
the exponent characterizing the stretched-exponential 
decay in regimes (b) and (c) varies slowly over time, 
due to crossover effects and amplitude corrections (fig- 
ure 1). It would be interesting to compare our theoreti- 
cal results to the anomalous decay of photogenerated ex- 
citons observed in TMMC chains and MX compounds. 
We note that the theory contains no free parameters 
beyond the diffusion coefficient of the excitons, D, and 
the concentration of disorder (defects), p. We have also 
considered a long-tailed distribution of segment lengths, 
4>(L) ~ L 1 (7 > 2). In this case the concentration 
decay is c ~ £ _ (t _1 )/ 2 ; faster than in the infinite line. 

The above results were extended to systems in higher 
dimensions, where the reaction zones could be mod- 
eled by percolation theory (e.g., annihilation of triplet- 
excitons in naphtalene-hs embedded in porous matri- 
ces of napthalene-dg). The results are similar to those 
in one dimension, but depend crucially on whether the 
system is above, below, or at the critical threshold of 
percolation, and d s (the spectral dimension of the clus- 
ters) affects the actual value of the various decay expo- 
nents. Below criticality, we find an early time decay of 
c ~ t~ d =/ 2 ~ i~ - 66 , followed by a stretched-exponential 
decay, c ~ exp[(const)i x ], with a slowly varying expo- 
nent: from x = d s /2 ps 0.66 to x = dg' /(2 + d s ') s» 0.37 
(d s ' is the spectral dimension of lattice animals, and all 
numerical values cited pertain to the physical case of 
three-dimensional percolation) |2q] . At criticality there 
is a faster power-law decay, c ~ £~ r C' r - 1 H/ 2 ~ t~ - 79 , and 
above criticality there is a crossover from c ~ £- d =/ 2 ~ 
£-0.66 a -|- ear jy times to a reaction-limited decay, c ~ t -1 , 
thereafter. 

As a final observation, we note that the scaling of 
the concentration decay for diffusion-limited annihilation 
(A + A — > 0) is exactly the same as that of coalescence, 
and hence the results presented here, in spite of being 
grounded on the exactly solvable case of coalescence in 
a segment, are valid also for annihilation. Likewise, the 
concentration decay at all times, and in particular the 
long-time exponential decay of eq. ([j~7|), remains essen- 
tially the same when the particles are absorbed at the 
segment edges, rather than reflected. Thus, our conclu- 
sions are pretty independent of the fate of the particles 
as they hit the boundaries of the reaction domains. 



APPENDIX A: RESIDUAL PARTICLES 

Consider a finite domain of volume V where coales- 
cence or annihilation might take place. At the end of 
the process, when no further reactions are possible, there 
would be left either one or zero particles. For coales- 
cence, A + A — > A, the residual number of particles is 
one, unless the volume had been empty to begin with. 
In the case of annihilation, A + A — > 0, the parity of the 
initial number of particles determines the residue: zero 
when the initial number is even, and one when it is odd. 

Suppose that the initial distribution of particles is com- 
pletely random, at concentration cq. The probability that 
the domain V contains exactly n particles is given by the 
Poisson distribution: 



Pn 



(c V) n 



-c V 



(Al) 



Thus, for coalescence, the probability that one particle is 
left over is 



p c I °f = l- Po = l-e 



-c V 



while for annihilation it is 



annil 
"res 



= Pn= 



1 - e- 2coV 



(A2) 



(A3) 



n odd 



If the system consists of disconnected domains of re- 
action, characterized by a distribution tp(V) of their vol- 
umes, and the initial concentration of particles is cq, then 
the concentration of residual particles at the end of the 
process is 



Cres = (V)- 1 / Ples (c ,V)lP(V)dV , 



(A4) 



where (V) — J Vip(V) dV is the average volume of the 
reaction domains. For example, if the distribution of 
reaction-domain sizes is exponential, ip(V) 



obtain 



„coal 



pe pV , we 



= cop/ (c + p) and 



c p/(2c + p). 



De Albuquerque and Lyra fll4| performed a lattice com- 
putation for the residual concentration of particles for 
annihilation in two-dimensional percolation, under the 
assumption that there is initially one particle per site. 
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